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(N| ■ ABSTRACT 

CO _ 

f*"^ ' This is the first of a series of papers in which the kinematics of disk galaxies over a range of scales is scrutinised employing spectroscopy. 

. A fundamental aspect of these studies is presented here: the new publicly available software tool TiRiFiC 
(http://www.astro.uni-bonn.de/~gjozsa/tirific.html) enables a direct fit of a "tilted-ring model" to spectroscopic data cubes. The algo- 
rithm generates model data cubes from the tilted-ring parametrisation of a rotating disk, which are automatically adjusted to reach an optimum 
Q ■ fit via a chi-squared minimisation method to an observed data cube. The structure of the new software, the shortcomings of the previously 
available programs to produce a tilted-ring model, and the performance of TiRiFiC are discussed. 

Our method is less affected by the well-known problem of beam smearing that occurs when fitting to the velocity field. Since with our method 
we fit many data points in a data cube simultaneously, TiRiFiC is sensitive to very faint structures and can hence be used to derive tilted-ring 
models significantly extending in radius beyond those derived from a velocity field. The software is able to parametrise H I disks of galaxies 
that are intersected by the line-of-sight twice or more, i.e. if the disks are heavily warped, and/or with a significant shift of the projected centre 
■ of rotation, and/or if seen edge-on. Furthermore, our method delivers the surface-brightness profile of the examined galaxy in addition to the 
k>( orientational parameters and the rotation curve. 
j_j . In order to derive kinematic and morphological models of disk galaxies, especially reliable rotation curves, a direct-fit method as implemented 
d ' in our code should be the tool of choice. 
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1. Introduction 

Kinematic analyses based on spectroscopy are an important 
tool to constrain the dynamical structure and hence the distri- 
bution of matter in galaxies. The discovery and description of 
global features in the k inematics of galaxies, such as the flat- 
ness of rotation curves (Rubin & Ford 1970; Bosma 1978), di- 
rectly influenced cosmology. Since then, the increase of com- 
putational power enabled theorists to provide testable predic- 
tions on the mass distributi on on sub-galaxy scales for given 
cosmological models ( e.g. iNavarro et alj 1997 ; Moore et al.l 
19991: iReed et al.l 120051) such that kinematical studies can be 
utilised as immediate tests for cosmology. While in recent years 
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the debate concentrated mainly on the sp herical distribution 
of Dark Matter in r elaxed systems (e.g . ISwaters et al.l I2003 



de Blok etal]|2003l: iGentile et all 12004 Navarro et alj \2004 



deviations from this, whic h are evident from lopsidedness (e.g. 
Schoenmakers et al.l l 1997) or warping (e.g. iBosmal 1978b . now 
come int o the focus of theoretical research in the cosmolo gical 
context dSharma & Steinmetzll2005l iGao & White] 120061) . The 
aim to constrain the structure of anisotropies in the DM distri- 
bution imposes an observational challenge, requiring both ob- 
servations with high sensitivity as well as appropriate analysis 
tools. 



This is the first paper in a series with the aim to push for- 
ward the observational limits in order to put further constraints 
on theory, to gain insight in the large-scale structure of the disk- 
halo system in disk galaxies, and to further test proposed mass- 
density profiles. One of the major ingredients is the necessary 
development of improved analysis methods. A new software 
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tool to analyse the kinematics of disk galaxies is presented in 
this paper. 

It is known from observations that the orbits of the disk ma- 
terial in spiral gala xies without large bars have a comparabl y 
low ellipticity (e.g. lBosmalll98ll ; ISchoenmakers et all 19971) . 
such that it is a good approximation to treat them as being cir- 
cular. This means that to first order the kinematics of a galactic 
disk at a certain galactocentric radius can be described by a 
set of three parameters, viz. the rotation velocity, and two pa- 
rameters that describe the local disk orientation with respect 
to some reference syste m. Such a "tilted-ring model" was first 
constructed for M83 by lRogstad et all dl974l) . 

In some cases the tilted-ring model is an oversimplifi- 
cation, for example when the orbits are sig nificantly non- 
circular due to th e presence of a bar (e.g. iBosmal 1 1978 



Simo n et al. I l2003l) . Nevertheless, in many cases it serves as 
a good approximation and hence is the standard kinematic 
model for galaxies. Several algorithms exist to fit such a 
model to spectrosc opic data. The most extensively used is th e 



ROTCUR routine dvan Albada etal.l |1985; Begemann 1987) 



which is a generalisation of the original method of lWarner et alj 
(119731) . It is implemented in several data anal ysis packages , 
e.g. GIPSY (Ivan der Hulst et ali i 19921). NEMO dTeubenll 19951) 
and MPS dFomalonllll98ll) . ROTCUR fits a set of inclined 
rings to a velocity field. ROTCUR derivatives and exten- 
sions exist. The GIPSY routine RESWRJ takes into account 
that non-axisymmetric potentials have a characteristi c imprint 
on the velocity field (iBinneyl Il978t iTeubenl Il991l) and al- 
low sjojjm^arimuj^l vamd (see 
also Franx et al. 1994; Schoenm akers et alj 119971) . Similarly, 
R1NGFIT dSimon et al.l 2003 ) allows for the same variation 



while at the same time only treating a flat disk. Si mon et al 



(2005) used this routine to test the reliability of their rotation 
curves. ROTCURSHAPE (implemented in NEMO, iTeuben 



1995) performs a global fit of the velocity field, in which the 
ring parameters are not fitted independently, but "in one go". A 
reduction of the higher order terms in a harmo nic expansion is 
used by KINEMETRY dKrajnovic et al.ll2006l) in order to find 
the best fitting solution of a tilted-ring model. 

All these routines have in common that they are based on 
an analysis of the velocity field, which itself is derived from 
a data cube, so it is an intermediate step when going from the 
observed data cube to a tilted ring model. Various methods ex- 
ist to extract a velocity field by analysing the spectra in a data 
cube. All of these methods have their shortcomings. First, for 
galaxies with large warps or galaxies seen close to edge-on, 
the derivation of a single representative velocity field is im- 
possible, because even if the disk geometry is known, the line- 
of-sight intersects the disk twice or more. This means that for 
some positions on the sky more than one velocity has to be in- 
ferred. Second, a velocity fie ld is contamin ated by beam smear- 
ing effects (see below, also ITeuben 2002 ). The integration of 
emission along the line-of-sight in a thickened galaxy causes 
a similar effect. This leads to a situation where, in principle, a 
tilted-ring model resulting from a fit to a velocity field has to be 
cross-checked with a model d ata cube as can be produc ed by 
the GIPSY routine GALMOD dvan der Hulst et all 19921 orig- 
inally designed by T.S. van Albada). Some authors have fit- 



ted such model data cubes by adjusting the model parameters 
and successively improving the model d ata cube by comparing 
it with the original data cube "by eye" dArnaboldi & Gallettal 
ll993HSwatersftl999tlGentile et al.ll2004h . 

The drawbacks of just using the velocity field fits led 
to the development of a number of software tools that per- 
form direct fits to the data cube (Irwin & Seaquist 1991) 
or to a position-velocity diagram dSimard & Pritchetl [l999; 



Takamiva & S ofue 2002; Bo hm et al.ll2004l) . thus circumvent- 
ing the beam-smearing or smoothing problem. These fit rou- 
tines, which simulate an observation from a model parametri- 
sation, do not allow, however, for an intrinsic change of the 
orientational parameters and hence are only limited implemen- 
tations of the tilted-ring model. Especially for the analysis of 
HI kinematics, full 3D fitting is most suitable, a nd has first 
been implemented by lCorbelli & Schneider! (119971) and applied 
to the spiral galaxy M33. In this case, however, the fit was per- 
formed on single spectra resulting from a single-dish observa- 
tion. Thus, besides not being available publicly, this algorithm 
would not suit to perform a fit to a data cube produced by syn- 
thesis or integral-field spectroscopic observations. 

In this paper, the new, publicly available software called 
TiRiFiC ('T/lted-TJ/ng-f/tting-Code") is presented, which per- 
forms an automated fit of a tilted-ring model to a data cube. 
This software was originally developed to construct tilted-ring 
models of heavily warped galaxies observed in the H I emis- 
sion line. It can, however, be used for all kinds of spectroscopic 
data cubes of translucent objects that can be approximated by 
the tilted-ring model. While the software will be under devel- 
opment for quite some time in order to improve the perfor- 
mance and to extend functionality beyond the classic tilted-ring 
model, TiRiFiC is well usable in its current form. 

In this paper, we try to answer two basic issues: does a 
method to fit a tilted-ring model directly to the data cube lead 
to more reliable results than a fit to the velocity field and is it 
thus possible to reach a reasonable fit within a reasonable com- 
putational time? The paper is laid out as follows: in Sect.|2]the 
model layout and the fitting procedure are introduced, in Sect. [3] 
the smoothing effect of velocity field-based fitting algorithms 
is discussed, in Sect. |4]the results of a few tests of TiRiFiC are 
presented. In Sect. [5] the results are summarised. First applica- 
tions of TiRiFiC to purposely selected (warped) galaxies will 
be published in forthcoming papers. 

2. TiRiFiC layout 

2.1. Tilted-ring model 

A TiRiFiC tilted-ring model is specified by a set of parame- 
ters that vary with radius, plus a set of global parameters. We 
refer to the set of parameters belonging to a given radius as a 
"ring" in order to stick to the traditional terminology. A model 
is calculated in a very similar, but not identical, way as in the 
GIPSY routine GALMOD. To calculate a model, a number of 
"sub-rings" with a user-specified width is created by linear in- 
terpolation of the ring-specific parameters. These sub-rings are 
then modelled by a Monte-Carlo integration to rotate with the 
same tangential velocity with a certain orientation w.r.t. the ob- 
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Fig. 1. Definition of the orientational parameters of TiRiFiC, viz. inclination i and position angle pa. Left: view of a circular orbit 
of the tracer material at an arbitrary position. A circular orbit appears to the observer as an ellipse as shown in the right graph. 
The three-dimensional orientation of the orbit circle is parametrised by the position angle enclosed by the north axis and the 
"descending" line-of-nodes (LON,desc), and the inclination. The descending line-of-nodes is the half-minor axis of the projected 
ellipse defined by the centre of the ellipse and the point where the galaxy material switches from having a higher recession 
velocity (than systemic) to a lower one, moving anti-clockwise along the ellipse. With the inclination being the angle between 
the celestial plane and the orbital plane, the orientation of the orbit is fixed. LOS is the line-of-sight, MA,app is the approaching 
half-major axis, MA,rec is the receding half-major axis, and LON,asc is the ascending line-of-nodes. 



server, determined by the position angle and inclination. The 
resulting velocities are then projected onto a cube with dimen- 
sions set by the input data cube. The orientational parametrisa- 
tion used by TiRiFiC is shown in Fig.Q] The final step to obtain 
a model representing an observation consists of a convolution 
with a 3D-Gaussian, representing the instrumental effects of a 
finite observing beam and the finite resolution in the frequency 
domain, as well as the internal dispersion within the tracer ma- 
terial of a galaxy. 

The user specifies the number of rings and for each ring 
(see also Fig.[T]): 

- the radius. 

- the circular velocity. 

- the scale height (for the specification of the vertical density 
distribution see below). 

- the surface brightness. 

- the inclination. 

- the position angle. Contrary to the usual convention (e.g. 
GALMOD, ROTCUR), the position angle is defined as the 
angle between the meridian and the minor axis of each 
ring from north through east (see Fig. [TJ. In comparison 
to GALMOD the TiRiFiC position angle is 90° larger . 
TiRiFiC uses the same definition as used by Briggs ( 1990). 

- the right ascension of the central positions of the rings, 
which is allowed to vary from ring to ring. 

- the declination of the central positions of the rings, which 
is allowed to vary from ring to ring. 

- the systemic velocity, which is also allowed to vary from 
ring to ring. 



The user specifies the following global parameters: 

- a global isotropic velocity dispersion, which includes the 
instrumental dispersion. It is clear that such a treatment 
is not entirely physical, as it is known that the dispersion 
changes with radius in spiral galaxies and is also not the 
same in the vertical and horizontal direction of the disk. 
On the other hand, one has to keep in mind that the tilted- 
ring model itself is in fact a geometric parametrisation of 
the galaxy kinematics, which still needs to be interpreted, 
even if it resembles the true kinematics quite well. The rea- 
son to keep a velocity dispersion as a global parameter is to 
improve the computational speed. It provides the possibil- 
ity to shift the modelling of the dispersion to a convolution 
instead of a (time-expensive) Monte-Carlo integration. 

- the second global parameter accounts for the vertical dis- 
tribution of the gas density. The user can chose between a 
Gaussian, sech 2 , exponential, Lorentzian and a box layer. 

- the constant total flux of a single point source (a "cloud"). 
The number of Monte-Carlo point sources is approximately 
the total flux of the model divided by this number. When 
calculating the point sources for a subring, however, the 
cloud flux is changed by a tiny fraction in order to keep 
the flux of the subring accurate. 

After calculation of a point source model and gridding of 
the Monte-Carlo point sources onto a model cube, the cube is 
convolved with a 3D-Gaussian, a product of a 2D (anisotropic) 
Gaussian in the x-y-plane, and a ID Gaussian determined by 
the global velocity dispersion. The form of the Gaussian in 
the jt-y-plane is determined by the observing beam of the in- 
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put cube (the CLEAN beam), which can be redefined by the 



user. 



2.2. x 1 evaluation 

The convolution routine implemented in TiRiFiC aims at 
computational speed. The only possibility to reach the re- 
quired speed is by an F FT-convolution, for which the FFTW 
( Frigo & Johnsonl2005l) library is used. It seems to possess the 
highest flexibility and quality (in terms of computational speed) 
of all freely available FFT libraries. 

TiRiFiC calculates the x 2 and, along with this, the relative 
probability of two models. Usually,^ 2 is calculated via 



X 



Z 



(M k - o k y 



-z 



(M k - Okf 



(1) 



where k is an index running over all pixels in the cube, M is 
the model data cube, O the original, and cr k the noise of the 
original data cube in the kth pixel. If the quantisation noise is 
treated correctly, the weight w k of a pixel should be 



w(k) = cr r 2 ms + (al) 2 , 



(2) 



where <T ms is the rms noise of the original data cube and 
cr q k the quantisation noise of the convolved artificial data cube. 
TiRi FiC is able to calculate cr^ (for a mathematical treatment 
J6zsall2006h . , however, a larger computational effort to 



see 



do so, and the user might be inclined to modify the goodness- 
of-fit evaluation. With a "weight parameter" W a weight map 
is calculated via 



w k = 



cr 2 • W 2 + (al) 2 



W 2 



(3) 



If W — > oo, the weight map becomes a constant, cr 2 ms , giv- 
ing TiRiFiC the possibility to save one convolution. If W = 1, 
the noise will be estimated correctly, including the quantisa- 
tion noise (which can in principle exceed the rms noise of the 
data cube). The parameter is kept continuous to give the user 
kind of a weighting scheme at hand. With an increasing weight 
parameter, the emphasis will be more and more taken away 
from regions with high quantisation noise (and thus high sur- 
face density), towards regions of low surface density. 

2.3. x 2 minimisation 

With the type of model generation in TiRiFiC, the usual x 2 
minimisation algorithms employing partial derivatives of the 
model function with respect to t he fitting parameters like the 
Levenberg-Marquardt a lgorithm dLevenbergl 1944 ; Marquardl 
1963c IPress et al.lll986f) are not easily realised, since the ana- 
lytic form of the fitting function is unknown. The Monte-Carlo 
method was chosen as a way to produce a model cube be- 
cause an analytic form of the tilted-ring model does not ex- 
ist. Furthermore, the model itself, for which a Jacobian would 
have to be supplied, is the convolved tilted-ring model in three 
dimensions. 

For TiRiFiC the simplest solution without the need for 
derivatives was implemented first hand. The golden-section 



search algorithm ([Press et al. 1986, Chap. 10.1) is an uneco- 
nomic but simple method (see below) to find local ^-minima 
in parameter space without the necessity to normalise the x 2 
evaluated with TiRiFiC, i.e. for this algorithm the knowledge 
of the rms noise in the analysed data cube does not play a ma- 
jor role. 

A tentatively impl emented Metr opohs Monte-Carlo- 



Markov-chain algorithm (iMetropolis et al.l l 1953") proved not to 
be applicable in tests using real H I data cubes, mainly because 
a smooth disk as fitted to the data cubes is not close enough to 
a realistic, clumpy galaxy disk. The Metropolis algorithm re- 
quires an ideal parameter-set to be able to reproduce the data 
perfec tly. For more details on the fitting procedure see Jozsal 
(120061 Chap. 2). 

Parameter fitting in TiRiFiC is rather flexible. The user can 
specify single parameters to be fitted, or groups of parameters, 
which will be fitted at once as one parameter. As an example, 
the user has the possibility to fit a changing projected centre 
from ring to ring as is, e.g., expected for a galaxy disk with a 
bowl-shaped warp. Such bowl-shaped warps result in an asym- 
metric projected appearance. For rather symmetric galaxies this 
might not be desirable, because in this case one would desire to 
keep the number of free parameters low. With TiRiFiC the user 
has the possibility to fit a single centre for all rings, defining 
their centres as a group and fitting the value of this group as a 
whole. 

3. Beam smearing and smoothing inherent to 
velocity fields 

The success of a tilted-ring fit to a velocity field depends on 
how well the velocity field is able to represent the mid-plane 
recession velocity of the rotating disk. In practice, to derive 
a velocity field from a data cube one usually tries to deter- 
mine the peak positions of the velocity profiles (e.g. ISwaters 
Chap. 3). The simplest way to extract this peak veloc- 



1999 



ity from a noisy velocity profile is to compute the inte nsity 
weighted mean (the fir st moment, e.g. iRogstad & Shostak 



1971 



(M± 



; Sim on et a 
Beaeman 




Garcfa-Ruiz et al 



2005b, whil e a fit of full Gaussians 



Swatersl 1_999) or half-Gaussians (e.g. 



2002) is considered a better approach. Even 
more sophisticated methods ex ist to reliably deriv e the peak 
velocities from the line profiles (Gent ileet al.ll2004 ). With arti- 
ficial data cubes, however, one is neither restricted by noise nor 
by the velocity resolution. Therefore, by producing a noiseless 
artificial observation with a very high spectral resolution and 
reading out the peak velocities from the profiles, one gets a ve- 
locity field which would ideally be derived from a measured 
data cube. 

It can be shown in a simple experiment that biases will oc- 
cur in a velocity field-based tilted-ring analysis even when us- 
ing such an idealised case of a velocity field. 

A Gaussian beam with a HPBW of 12 " and a (total) veloc- 
ity dispersion of lOkms -1 was used to produce with TiRiFiC 
two low-noise artificial observations of a galaxy with 17.6 • 10 6 
point sources in a data cube with a spatial pixel separation of 
4" (for the parametrisati on see Fig.|2l a detailed description can 



be found in Jozsa 2006). The parameters were sampled with a 
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arcsec 



Fig. 2. Testing the extraction of rotation curves from velocity fields. Filled circles and connecting lines (black) show the parametri- 
sation with which two artificial observations have been generated, a velocity field and a data cube (see text). Top panel: rotation 
curve. Centre panel: blow-up of the rotation curve (same lay-out as above). Bottom panel: H I surface brightness. The inclinations 
are 60° (left) and 75° (right) respectively. Velocity fields were generated as described in the text. Triangles and connecting lines 
(blue in online version): results from a ROTCUR fit with a free-angle of 0° and a uniform weighting. Stars and connecting lines 
(red in online version): results from a ROTCUR fit with a free angle of 30° and a cosine weighting. Diamonds and connecting 
lines in the right panel (green in online-version): results from a cut along the galaxy major axis, corrected for inclination. Filled 
circles and connecting lines (black) also denote the results of a TiRiFiC fit; they are identical to the parametrisations of the 
artificial observation. 



separation of 12". Inclinations of 60° and 75° were chosen, and 
the position angle was kept fixed at all radii. 

In order to produce a "perfect" velocity field, the cubes 
were generated with an unrealistic channel separation of 
0.26 km s~'. With these data cubes peak-velocity maps were 
created. Using these velocity fields, rotation curves were pro- 
duced using the GIPSY routine ROTCUR, employing two 
slightly different approaches with an input guess identical to 
the model parameters, using a sampling of the rings identical 
to the inp ut model. The results of both fits can be seen in Fig. [2] 
(see also lTeubenll2002l) . 



The fitted rotation velocities are biased following a certain 
pattern that depends on the shape of the surface brightness pro- 
file, which in this test case is falling monotonically. In those re- 



gions where the rotation curve is rising, the rotation velocities 
are underestimated, while in regions where the rotation curve 
is falling the rotation velocities are overestimated. Hence, the 
fitted rotation velocity tends to be biased towards regions of 
higher surface brightness. The reason for this lies in the fact 
that the measured velocity profile at a certain position of the 
galaxy is a convolution of the beam with the true intensity dis- 
tribution. Thus, if the intensity is not constant over the beam 
area, the resulting "smeared" profile receives more contribu- 
tion from areas with higher intensity and hence becomes asym- 
metric, even if it would be symmetric when observed with a 
pencil-beam. The position of its maximum is dragged towards 
velocities with a higher associated surface brightness. This ef- 
fect is cannot be treated in a straightforward manner because it 
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also depends on how velocities change within a beam. It is also 
visible at radii where both the rotation curve and surface bright- 
ness are flat over a large range. This means that "smoothing" 
affects the rotation curve at all radii. The reason for this lies in 
the fact that in projection the galaxy becomes fore shortened 
along the kinematical minor axis. Here, the gradient of the ve- 
locity field becomes largest and the observing beam covers the 
largest range of velocities along the kinematical minor axis. 
Points in the vicinity of the kinematical minor axis are thus 
more strongly affected by beam smearing. This is why usually 
a section of the velocity field is thus not regarded (by specify- 
ing a "free-ang ') in a rotation curve analysis and a weight- 
ing of data points is applied, while our experiment shows that 
this treatment does not help much for the analysis of the inner 
regions, but improves matters in the outer ones. An extreme 
choice for a free angle would be to only read out the velocities 
along the major axis of the galaxy and to correct the velocities 
for the inclination. This is equivalent to an extraction of a rota- 
tion curve from a long-slit observation. The results from such 
a "fit" are shown in the right hand panel of Fig. [2] It can be 
seen that with such a choice the results get slightly better in the 
outskirts of the galaxy at the expense of far worse results in its 
inner parts. 

In order to demonstrate how a velocity field analysis is in- 
fluenced by the a-priori unknown distribution of the Hi, we 
chose a rather high inclination of at least i = 60° for our artifi- 
cial observations. Since this covers only one third of all observ- 
able galaxies, we emphasise that not necessarily all velocity 
field analyses are affected by beam smearing. Furthermore, the 
dumpiness of the distribution of the interstellar medium (ISM) 
introduces further errors in a tilted-ring analysis, which might 
well exceed the ones arising from the shown variant of beam 
smearing. 

However, using a velocity field as would be ideally derived 
in a noiseless observation, we demonstrated how in analysing 
velocity fields of galaxies with an inclination higher than i = 
60° systematic errors might be introduced by beam smearing. 
Even if the maximum deviation from the true rotation curve 
is not very large (about 3.3kms _I for i = 60° and about 
3.9kms~' for i = 15° in this experiment), this variant of beam 
smearing can introduce a significant misinterpretation of the 
measurement, as the resulting rotation curve is systematically 
biased. Moreover, second order treatments of the kinematics of 
gaseous disks, like a harmonic analysis, require a precision in 
that range. An assignment of a (large) statistical error to com- 
pensate for beam smearing is thus not a solution. Results from 
mass decompositions in which (unknown) density profiles are 
fitted to rotation curves derived from velocity fields are very 
vulnerable in that respect. The smoothing effect can lead to a 
preference of rotation curves with a smaller curvature d 2 v/dr 2 . 



1 angle defining cones around the minor axis within which radial 
velocities are discarded from the computation 



4. Testing TiRiFiC 

4.1. The trivial test 

Figure [2] also shows the results of a fit with TiRiFiC to arti- 
ficial data cubes that are identical to the data cubes used to 
produce the velocity fields, except for a larger channel width of 
2.06 km s~'. As for ROTCUR, the first-guess-parameters were 
identical to the model parameters and sampled at identical 
radii. Naturally, TiRiFiC makes a perfect fit and the resulting 
model is identical to the input model. It is however noticeable 
that, in principle, with TiRiFiC it is thus possible to produce an 
artificial observation that can be reproduced by a fitting process 
over and over. For routines working on the velocity field this is 
not the case. First, the surface brightness distribution, which is 
a free extra output of TiRiFiC, has to be inferred using other 
information (e.g. by fitting a tilted-ring model to the total in- 
tensity map). Second, an artificial data cube will always have to 
be subject to a convolution with an artificial beam, hence intro- 
ducing beam smearing; therefore, the artificial cube never re- 
produces a velocity field that fits to the tilted-ring model found 
by the fitting routine. 

4.2. Rotation curve of an edge-on disk 

One of the strengths of TiRiFiC is supposed to lie in the abil- 
ity to perform reliable fits where velocity field methods are 
bound to fail by definition. This is in particular the case when 
the line-of-sight crosses the galaxy disk twice or more often. 
An extreme case of such a situation is encountered when a 
galaxy is observed (nearly) edge-on, because then the galaxy 
disk is crossed not only twice but an infinite number of times. 
We simulated a galaxy observation using the rotation curve and 
the surface brightness profile shown in Fig. [3] a template cube 
and an artificial beam as in Sect. I4.ll putting it in an edge-on 
orientation. The parameters for the artificial observation were 
sampled in steps of 6". Then, the artificial observation was fit- 
ted with TiRiFiC, taking the input model parameters as a first 
guess, but leaving the rotation velocity, the surface density pro- 
file, and the orientational parameters, inclination and position 
angle as free parameters. In order to prevent a trivial solution, 
the fitting was performed with a parameter sampling in steps 
of 12". TiRiFiC reaches a nearly perfect fit, interpolating in 
the regions where the fit results cannot reproduce the artificial 
observation by construction (see Fig. [3). Hence, TiRiFiC is in 
principle able to figure out reliable rotation curves from obser- 
vations of edge-on galaxies. A slight "beam-smearing" effect 
takes place at radii where it is not possible to follow the high 
curvature of the parametrisation of the artificial observation, 
owing to the fact that in the fitting process the parameters are 
sampled with half the rate as for the artificial observation. 

A caveat is that TiRiFiC assumes a galaxy to be translucent. 
Observations of edge-on galaxies may suffer significantly from 
self- absorption and scattering processes at certain wavelengths 
(e.g. Bae s et al.ll2003l) . This has to be taken into account when 
interpreting TiRiFiC results, especially when fitting to optical 
spectroscopic data. 
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Fig. 3. Testing rotation curve extraction with TiRiFiC. 
Triangles and connecting lines (red in online version) show 
the parametrisation with which a artificial observation of a flat 
galaxy with 90° inclination has been generated (edge-on case). 
Top panel: rotation velocity. Middle panel: rotation velocity, 
blow-up. Bottom panel: surface brightness. The parametrisa- 
tion for the artificial observation is twice as dense as for the fit 
of the artificial observation performed with TiRiFiC. The fitted 
parameters are represented by dots (black in online version) 
and connecting lines. 



Fig. 4. Testing tilted-ring fitting with TiRiFiC. Triangles and 
connecting lines (red in online version) show the parametri- 
sation with which a artificial observation of a heavily warped 
galaxy has been generated, stars and connecting lines (blue in 
online version) show the first-guess input to TiRiFiC. The fitted 
parameters are represented by black filled circles and connect- 
ing lines. 



4.3. An edge-on warp 

TiRiFiC has been constructed to be able to fit data cubes of 
observations of warped galaxies, which is one of the main ap- 
plications of the tilted-ring model. In order to construct a case 
tailored to put TiRiFiC to the test and for which no appropri- 
ate automated software exists yet, a warp in an edge-on galaxy 
was constructed. This is a situation where it is impossible to 
construct an unambiguous velocity field, as the line-of-sight 
crosses the disk several times. Again, the model observation 
was generated using a template data cube and an artificial beam 
as described in the previous subsections, sampling the model 
parameters in steps of 6". 

The input parameters to TiRiFiC were chosen to randomly 
deviate from the input model parameters (see Fig.|4|i. Again, in 



fitting the data the model parameters were sampled in steps of 
12" in order to prevent the possibility of a trivial solution. 

In fitting the data cube, the rotation velocity, surface bright- 
ness, inclination and position angle were fitted independently 
for every single ring; the other parameters were fitted by letting 
them vary together. An exception was made for the innermost 
four rings where the orientation parameters were fitted with the 
constraint to be the same for each ring. A number of 3.5 ■ 10 6 
point-sources was chosen for the fitting process, while the arti- 
ficial observation was created with 17.6 • 10 6 point sources. 

The result of the fit process (Fig. @]l shows that with 
TiRiFiC, it is possible to fit to a data cube of an edge-on galaxy 
with a large warp with a very high accuracy. 
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4.4. Sensitivity 

Since a galactic disk (or any gaseous astronomical object with 
disk symmetry) becomes fainter with increasing radius, the ex- 
tent of a tilted-ring analysis is either limited by the observa- 
tional noise or by the fact that the ring symmetry breaks down. 
For a conventional velocity field analysis the radial velocities 
have to be derived from the spectra in the single pixels in a 
data cube. Hence, the extent of the velocity field and, with this, 
the maximal radius that can be addressed in a tilted-ring anal- 
ysis is determined by the pixels where the signal-to-noise ratio 
significantly exceeds unity. Since with the TiRiFiC approach 
one fits a global model to the data cube, the number of pixels 
that determine the parametrisation of a ring is large. Therefore, 
even if the emission from an object is well below the nomi- 
nal noise per pixel, one should still be able to derive a reliable 
parametrisation with TiRiFiC . Provided the object follows a 
tilted-ring symmetry the detection limit should in theory scale 
with 1 / V(A0, where N is the number of pixels with emission 
in a rotating ring projected onto the 3D data cube. 

In order to get an impression of the realistic performance 
of TiRiFiC when fitting faint structures, we generated a large 
number of artificial observations of gaseous rings with vary- 
ing uniform surface brightness and varying surface solid an- 
gle. The rings were projected onto a data cube with a size of 
512" x 512" and a pixel size of 4", 64 channels with a separa- 
tion of 4.12kms~' using an observational beam of 12" x 14" 
(Half Power Beam Width), and a global velocity dispersion of 
7 km s . With these specifications, an H i observation with the 
Westerbork Radio Synthesis Telescope (WSRT) is simulated. 
In order to add realistic noise to the model, we utilised line-free 
channels from a 2 x 12 hrs observation with the WSRT having 
an rms noise level of 3.6 mJy/beam. The rings were chosen to 
have an inclination of 60°, a rotation velocity of 75 km s , and 
a scale height of 2" using a sech 2 -law for the vertical distribu- 
tion. While the maximal extent of each ring was fixed at 200", 
the inner radius was varied between 45" and 195", varying lin- 
early the face-on solid angle occupied by the ring, which is 
roughly proportional to the number of pixels containing emis- 
sion in the projection of the ring onto the cube. 

We fitted the artificial observations with rings of the 
same extent, with an initial guess deviating significantly 
from the artificial parametrisation (surface brightness 1 • 
10~ 5 Jykms -1 arcsec~2, rotation velocity 65kms~', inclina- 
tion and position angle deviating by 10°, all other input pa- 
rameters identical to the parametrisation of the artificial ob- 
servation) in a two step process, first fitting all parameters ex- 
cept the surface brightness, then fitting all parameters. In order 
to distinguish a detection of such a ring from a non-detection 
we also fitted a data cube containing only noise and no emis- 
sion. Calculating the differences of the fitted parameters and the 
parameters used for the artificial observation, one can clearly 
distinguish between cases where the fit was successful and 
where it failed (Fig. [5] and Figs. [8] in the online version of 
this paper). By setting a threshold for the deviations of the pa- 
rameters, viz. rotation velocity (4kms _1 ), surface brightness 
(7- 10~ 7 Jy km s _1 arcsec~2), position angle (2°), and inclination 
(2°), a binary diagram could be constructed in which a detec- 




2e4 4e4 6e4 8e4 10e4 12e4 

ring solid angle (sq.arcsec) 

Fig. 5. Testing the sensitivity of TiRiFiC, example. Each pixel 
represents a fit to a artificial observation (including realistic 
noise) of an inclined rotating ring with a given solid angle and 
surface brightness, the lowest row representing a fit to a data 
cube containing noise only. The pixel values represent the de- 
viation of the fitted position angle from the true position angle 
of the ring. Clearly distinguishable are regions where the fit 
succeeded (bright regions) and where it failed (dark regions). 
Similar plots can be generated for other fitting parameters (see 
Fig. [8] in the online version of the paper). 

tion for combinations of ring surface brightness and ring solid 
angle was marked (see lower right panel in Fig. [8j. Finally, the 
detection limits as a function of the ring solid angle could be 
extracted (Fig. |6j. It turned out that while the detection limit 
determined that way shows a monotonic decrease with increas- 
ing ring solid angle, the expected proportionality did not show 
up. Also when trying to estimate the number of independent 
data points which the projected ring occupies in the artificial 
observation we could not establish a simple law that connects 
the ring solid angle to the detection limits. 

However, Fig. [6] clearly shows how the sensitivity in de- 
tecting faint symmetric structures in gaseous disks is enhanced 
using the TiRiFiC approach. With a power law 



q-0.82 



(4) 



where cr t ; r is the detection limit and £2 r i ng is the ring surface 
solid angle, and £2 m b = 202 arcsec 2 the beam solid angle, we 
achieve a reasonable fit. With this we can estimate the ring solid 
angle of 141 £2 m b that is needed to detect a ring that in projec- 
tion onto the cube has an intensity just below the noise <x rms in 
the data cube (1.6 • lO^Jykms -1 (&2 ■ 10 19 atoms cirT 2 in an 
Hi observation, horizontal line in Fig.|6]l. At a radius of 200", 
this corresponds to a ring width of 24". 

A caveat arises from the fact that even when no real emis- 
sion is present, a fit with TiRiFiC will result in a detection a 
faint structures in the presence of noise: in a noisy data set a 
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Fig. 6. Testing the sensitivity of TiRiFiC. Detection limits as a 
function of ring-surface solid angle in units of the main beam 
solid angle (Q m b = 202arcsec 2 , see text). The dependence is 
fairly well fitted by a power law (curved line). The vertical line 
marks the surface brightness where all pixels in the noiseless 
data cube lie below the l-cr rms level of the noisy data cube. The 
right hand axis is scaled such that the surface brightness of a 
ring is converted to surface column density of an H i observa- 
tion. 



residual offset in surface brightness can be surely found for 
a specific optimal ring parametrisation even if averaging over 
many pixels. This can be tested when fitting to a data cube that 
contains only noise, like it has been done in the experiment dis- 
cussed here. Figure [7] shows the detected surface brightnesses 
as a function of ring width when fitting a data cube that con- 
tains only noise. The erroneously detected surface brightnesses 
lie well below the detection limit as quoted in Eq. 01 the de- 
tection limit exceedig the erraneously detected surface bright- 
nesses by a factor of a few. However in analysing observations 
with TiRiFiC, one has to be cautious when accepting a faint 
structure detected by TiRiFiC as real, since then residual cali- 
bration and CLEANing errors become more and more impor- 
tant; in particular cases one should possibly apply additional 
tests such as the analysis of independent observations. 

Compared to e.g. H i observations of nearby spiral galaxies, 
the size of the test data cube is rather small. Since the surface 
solid angle of a ring with constant width scales with radius, 
the detection limits will scale with radius for larger objects. 
We thus conclude that with TiRiFiC, one is able to detect faint 
structures well below the noise - provided the observed object 
shows the symmetry inherent to the tilted-ring model. The de- 
tection limits can be estimated as shown above. 

4.5. Realistic tests 

In order to test how TiRiFiC works in practice, we constructed 
a number of artificial H i observations of galaxies that had to 
be analysed with TiRiFiC by a (variable) author of this paper, 
the "observer". The template cube and the observing beam had 
the same specifications as given in Sect. 14.41 We performed fits 
to artificial observations with and without realistic noise added 
to the data cube as described above, and with varying degrees 
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Fig. 7. Resulting surface brightness from a fit of a ring to a cube 
that contains only noise as a function of ring surface solid angle 
in units of the beam solid angle (£2 m b = 202 arcsec 2 , see text). 
Over-plotted is the detection limit determined from fits to faint 
rings (curved line, see text). The erroneously detected surface 
brightness lies well below that line. 



of symmetry. The observer was informed about the symme- 
try of the object (like "flat disk" or "constant centre and sys- 
temic velocity"), but strictly not about the specific parametrisa- 
tion of the artificial galaxy disk. We also included cases where 
all parameters were kept variable and hence a lopsided or a 
U-shaped galaxy with a possible flare was analysed. Finally, 
we constructed a case where the cloud flux was drastically en- 
hanced in order to simulate the dumpiness of the gas in the 
observed object. 

All fits were performed following the same pattern. First, 
the observer estimated the rough geometry of the object from 
the data cube, in order to then run a first iteration process with 
TiRiFiC using only a few parameters (fitting a common cen- 
tre and systemic velocity, a common scale height, a common 
orientation, and a surface brightness distribution and a rota- 
tion curve parametrised by a few data points). After that, the 
sampling of the parameters was refined and, if necessary, the 
parameters were fitted independently. It turns out that a fitting 
process with TiRiFiC needs the personal attention of an ob- 
server. In order to achieve a good fit, it is necessary to inspect 
the results and to correct single data points that are obviously 
"outliers" in order to restart the fitting process. After a few such 
iterations (in the range of 2-10) no improvement is achieved 
and the last fit is the final result. 

In order to compare TiRiFiC to methods, where a velocity 
field is derived and fitted, i.e. a ROTCUR application together 
with an idealised realistic construction of a velocity field, we 
additionally generated data cubes with a ten times higher ve- 
locity resolution (a channel width of 0.41 kms~') and (if the 
cube given to the observer contained noise) a ten times lower 
noise level. The data cubes were flagged where the intensity in 
the noiseless data cube was below 0.5<x rms . Both a peak inten- 
sity and a first moment velocity field was generated from the 
flagged data cube using the GIPSY routine MOMENTS. This 
way we simulated the construction of a velocity field as would 
ideally be derived by an observer (see Sect. 0) using common 
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methods for an extraction of a velocity field. The velocity fields 
were analysed with ROTCUR, the input guesses either deviat- 
ing slightly from the optimal guess or being optimal. 

A full documentation of all specific fitting results is given 
in the online version of the paper. Here, we provide a summary. 

Independent of the specific iteration method preferred by a 
single observer, TiRiFiC bears reliable results at galactocentric 
radii greater than 1-1.5 HPBWs (in our example greater than 
20")- Below that radius, the number of independent data points 
is too low to allow a fit of the full parameter set. In a statistical 
sense, TiRiFiC is thus not free from beam smearing but by far 
less affected than the traditional method: for the bulk of the data 
points, the ROTCUR results are worse than those of TiRiFiC, 
and the minimal radius where ROTCUR produces reliable re- 
sults (under idealised circumstances) lies beyond 4-5 HPBWs 
of the observing beam. 

From Eq. |4] we estimate the detection limit for a ring 
width of 15" at a radius of 200" to be 2.3 • 10~ 5 Jykms" 1 
(& 2.9 • 10 19 atoms cm for an H i observation). This value is a 
conservative upper limit. In most cases the detection limits are 
significantly lower, especially when the observer could make 
use of a certain symmetry, e.g. a common centre of all rings. 

TiRiFiC fits reliably disks that possess shifting centres and 
systemic velocities. Hence, with TiRiFiC one is able to fit 
asymmetric galaxies that maintain the tilted-ring symmetry. 

A clumpy Hi distribution simulated as described above 
does not hamper the performance of TiRiFiC. TiRiFiC aver- 
ages over azimuth, hence a random distribution of a still size- 
able number of clouds results in a larger^ 2 , but not in a differ- 
ent result. We did not simulate the case where the Hi distribu- 
tion is distinctly asymmetric. It is to be suspected that in some 
cases an m=2 harmonic in the distribution of the tracer material 
can have an effect on the results. 

Above a limit of 25° in inclination (which is close to (face- 
on) TiRiFiC disentangles reliably rotation velocity and inclina- 
tion. In this respect, ROTCUR made a better performance , due 
to to the fact that our derived velocity fi elds become very accu- 
rate at low inclination. It was shown by lBegemannld 19871) . that 
for an inclination lower than 60°, inclination and rotation veloc- 
ity become significantly degenerate in a tilted-ring analysis of 
the velocity field. For a velocity field analysis with ROTCUR, 
to derive reliable rotation curves (and with that, a reliable es- 
timate for the inclination) of galaxies with an inclination of 
less than 40° it is necessary to apply the method on a veloc- 
ity field with a statistical error of less than 1 rmkms' 1 . Such 
an accuracy is hard to reach using real observational material. 
In cases with low inclination, the surface density profile deter- 
mined with TiRiFiC had the smallest errors. 

For most parameters a stable solution is found in the first 
iteration process. The fact that nevertheless a few iterations are 
needed to reach a final model shows that for some singular pa- 
rameters the algorithm used for the x 2 minimisation does not 
perform an effective scann of the parameter space. Then, the 
user takes the role of stirring up the parameters in order to 
enable the localisation of the better minimum (for a few data 
points). The success and the quality of the fits in our tests per- 
formed was independent of identity of the observers, who used 
slightly differing approaches with regards to the fitting parame- 



ters and the choice how to approach the best-fit parametrisation 
from start on. Thus, it seems that the actual way the user shifts 
singular data points in order to start a next iteration process has 
no big influence on the quality of results. 

We fail to achieve good results in one particular case with 
both, TiRiFiC and ROTCUR. In the case of a galaxy with a 
thick disk (a vertical scale height of 33" in a disk with an Hi 
radius of about 200", see Fig.[l0]and Fig.[TT]in the online ver- 
sion of the paper), a TiRiFiC fit will converge extremely slowly, 
since in this case the degeneracy of inclination and rotation 
velocity is enhanced. This slow convergence then leads to a 
misinterpretation of the results as being already the final solu- 
tion. Also ROTCUR fails with extremely large errors, since the 
kinematical information is smeared out along the kinematical 
minor axis, making the x 1 minimum in parameter space shal- 
lower. In other words, the iso-velocity contours get stretched 
artificially along the minor axis. In case of TiRiFiC we have 
the hope to overcome this problem by implementing a more ef- 
fective minimisation method. The user will in any case be able 
to detect a thickened disk, and be aware of the problem. We ex- 
pressively caution, however, against the over-interpretation of 
rotation curves regardless of the analysis method in the case of 
galaxies likely to have a thickened disk like it is t he case for dlrr 



galaxies (see e.g. [S tavelev-S mith et al.l 11992c iBottema et al 



ll986hlWalter & Bri nks 1999). Without showing a galactic disk 
to be thin enough, these rotation curves are extremely unreli- 
able. 



4.6. Solid-body rotation 

Another problem occurring when extracting rotation curves us- 
ing velocity fields occurs in the case of a soli d-body rotation, 
as ca n be found again in dlrr galaxies (e.g. IWalter & Brinks! 
1999b . In that case, rotation velocity and inclination are com- 
pletely degenerate. Since with TiRiFiC one fits not only the 
kinematics but also the radial surface brightness profile, this 
degeneracy is broken, provided the gas distribution is symmet- 
ric enough and shows a radial gradient, as we can show in a 
further experiment (last viewgraph in the online version of the 
paper) that we performed with the same specifications for the 
artificial observation as given in Sect. l4.5l using a noiseless data 
cube. TiRiFiC fits very reliably a flat disk with solid-body ro- 
tation. 



4.7. TiRiFiC performance 

TiRiFiC is slow. Nearly 160 hours computation time needed to 
perform the fit described in Sect. 14.31 working with an AMD 
Athlon XP 2400+ CPU, can serve as an orientation. In prac- 
tice, however, the user has the possibility to reduce the com- 
puting time considerably. In a first iteration process, the size 
of the data cube can be reduced by binning data points. The 
model can be ported easily from a low-resolution data cube 
to a high-resolution data cube, since the model parameters are 
independent of the specifications of the data cube, unless the 
reference coordinate system (the epoch and the velocity defi- 
nition) changes. Assumptions about symmetries can be made 
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in order to reduce the number of parameters. E.g., if a disk is 
not warped, all position angles and inclinations can be fitted as 
one single parameter. A-priori known parameters like a surface 
brightness profile can be excluded from the fitting process. 

Up to now, we concentrated on an implementation of a 
working version of the software in order to demonstrate its ap- 
plicability. Future development, however, will include a paral- 
lelisation of the code in order to make the use of computer clus- 
ters possible. Furthermore, the possibilities to use more sensi- 
ble fitting procedures will be explored. With this we certainly 
expect a jump in the computing performance of TiRiFiC. 

5. Results and discussion 

We introduce TiRiFiC as a program that fits a tilted-ring model 
directly to spectroscopic data cubes and we provide a num- 
ber of tests of the program. We show that the TiRiFiC tilted- 
ring model has a high flexibility and even allows modelling of 
galaxy disks that are asymmetric in projection. We are able to 
show that TiRiFiC is not affected by the drawbacks of methods 
that work on velocity fields. For galaxies with an inclination 
above 60°, we show that velocity fields are affected by sys- 
tematics that depend on a-priori unknown factors, such as the 
surface density distribution. While the deviations of a rotation 
curve from the true one derived with fits to the velocity field are 
small, they occur in a systematic manner, possibly biasing the 
results of a rotation curve analysis. Using a direct fit to the data 
cube, TiRiFiC reduces these systematics significantly, deliver- 
ing the unknown surface density profile as an extra. A finite 
observing beam affects the reliability of the TiRiFiC results at 
the smallest radii in a statistical sense, since the number of fit 
parameters becomes too large in comparison to the number of 
independent data points. In general, TiRiFiC delivers more re- 
liable results than a fit to the velocity field. 

The biggest advantage of TiRiFiC lies in its applicability to 
edge-on and warped galaxies. Several tests show that TiRiFiC 
is able to reach a best-fit model that comes very close to artifi- 
cial input models having a large edge-on warp. 

Since the approach of TiRiFiC is to fit a global model to the 
data cube, a detection and quantification of a faint outer disk is 
possible, even if the bulk of the emission of that disk is fainter 
than the rms noise of the data cube, provided the outer disk 
follows the tilted-ring symmetry. 

TiRiFiC fails to disentangle inclination and rotation veloc- 
ity below a disk inclination of 25° (i.e. close to face-on). In that 
case, to estimate a more reliable rotation curve, a velocity field 
approach is to be preferred. 

The determination of rotation curves for galaxies with very 
thick disks like dlrr galaxies is unreliable, regardless of the fit- 
ting method, since the degeneracy of inclination and rotation 
velocity is significantly enhanced. Since the surface-brightness 
distribution is also fitted by TiRiFiC, it is nevertheless possi- 
ble with this approach to derive reliable rotation curves in the 
case of solid-body rotation that occurs predominantly in din- 
galaxies. 

The major current drawback of TiRiFiC is its slowness. 
Performing a single fit when using a conventional PC takes a 
long time. Hence, TiRiFiC is not (yet) usable to perform many 



fits in a short time unless a certain amount of computational 
power can be employed. 

We present a series of practical tests in the course of which 
an "observer" had to parametrise an unknown galactic disk 
with TiRiFiC, and we are able to show that such a parametri- 
sation is reliable even when analysing data cubes with real- 
istic noise. The success in the limits stated above is indepen- 
dent of the specific approach of the observer. It is of course 
the application to real observational data as will be presented 
in subsequent papers of this series that will show how useful 
TiRiFiC is in gaining new insights in the kinematic structure 
of rotating disks, especially galaxy disks. TiRiFiC can be con- 
sidered as the most straightforward realisation of a tilted-ring 
fitting method, since with its approach one does not require any 
careful analysis of the single spectra in a data cube in order to 
extract a velocity field. Hence, if TiRiFiC fails to produce a 
model data cube that fits to an observation of a galactic disk, 
the implication is that its kinematics cannot be represented by 
a tilted-ring model. The development of fitting routines similar 
to TiRiFiC is a good step towards improving analysis meth- 
ods of observational data containing kinematical information, 
as will be demonstrated in the consecutive papers of this se- 
ries. As secondary kinematical features (and their influence on 
the analysis) not compliant with the tilted-ring model have re- 
ceived more and more attention in the recent literature, e.g. by 
applying harmonical analyses, the aim must be to extend the 
analysis tools in order to put theory of galactic structure forma- 
tion that includes non-circular motions to the test. This work is 
already in progress. 
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6. TiRiFiC sensitivity: additional plots 
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Fig. 8. Testing the sensitivity of TiRiFiC. Each pixel represents a fit to a artificial observation (including realistic noise) of an 
inclined rotating ring with a given solid angle and surface brightness, the lowest row representing a fit to a data cube containing 
noise only. The pixel values represent the deviation of the fitted surface brightness (upper left panel), the rotation velocity (upper 
right panel), and the inclination (lower left panel) from the known true values. Clearly distinguishable are the regions where the 
fit succeeded (bright regions) and where the fit failed (dark regions). By setting a threshold a binary diagram can be generated 
(lower right panel), from which the detection limits in dependence of the ring surface solid angle can be estimated (Fig. [6]). 
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7. Testing TiRiFiC: individual tests captions. The fit was performed using a free angle of 45°, fix- 

ing the expansion velocity of the rings to Okms . For each 
test we provide a plot (left) containing the parametrisation of 
the artificial disk (black dots and connecting black lines), the 
results of the TiRiFiC fit (red triangles and connecting lines), 
and the results of both ROTCUR fits (blue diamonds and con- 
necting lines: fit to first-moment velocity field, stars in magenta 
and connecting lines: fit to the peak velocity field). From top 
to bottom the panels contain: SD: Hi face-on column density 
and surface density in atoms cm -2 and M pc~ 2 , respectively. 
VROT: rotation velocity in kms -1 . INCL: inclination in de- 
grees and radians, respectively. PA: position angle in degrees 
and radians. RA: the right ascension of the central position in 
hh:mm:ss. RASH: the projected shift in right ascension in kpc 
with respect to the centre of ring 1. DEC: the declination of 
the central position in dd:mm:ss. DESH: the projected shift in 
declination in kpc with respect to the centre of ring 1. VSYS: 
the systemic velocity in kms 1 . SCHT: the scale height in arc- 
seconds and pc, respectively. The black vertical line marks the 
end of the artificial disk, while the red vertical line marks the 
end of the disk modelled by the observer. The blue vertical line 
marks a radius of 1.5 HPBWs, beyond which TiRiFiC pro- 
duces reliable results in nearly all cases. The blue horizontal 
line marks the detection limit for a ring of 15" width at a radius 
of 200" as derived in Sect. 14.41 The right-hand graph shows 
the deviations of the TiRiFiC fits from the model parametrisa- 
tions. Black dots and lines represent the deviation of the black 
dots (the true parametrisation of the artificial galaxy) in the left 
viewgraph from the red lines (the parametrisation fitted with 
TiRiFiC) in the left panel, in red the deviation of the red trian- 
gles (the parametrisation fitted with TiRiFiC) in the left panel 
from the black lines (the true parametrisation of the artificial 
galaxy) in the left panel. Blue diamonds show the deviation 
of the ROTCUR fit to the first-moment velocity field from the 
parametrisation of the artificial galaxy, stars in magenta show 
the deviation of the ROTCUR fit to the peak-velocity field from 
the parametrisation of the artificial galaxy. The vertical lines in 
the right-hand plot have the same meaning as in the left-hand 
plot. 

We analysed the quality of the TiRiFiC fits by means of 
statistics. To this end, we extracted the deviations of the fits 
from the parametrisation of the artificial galaxy and built the 
mean and rms, and calculated the maximal error. The tables 
contain: 

(1) The quantity listed. A: mean of deviations. VP: rms of de- 
viations. A max : Maximum deviation in the range described by 

(2) . 

(2) The range in which the deviations were determined, all: 
taking all data points into account. 70%: taking the 70% low- 
est deviations into account, n - cr. taking those data points 
into account where the expected surface brightness is below 
the n - cr- level of the artificial observation. The 1 - cr level 
is determined by the face-on column density of a ring con- 
structed as described in Sect. 14.41 which in projection onto a 
data cube contains no emission above the 1 - cr rms level in the 
noisy data cube (cr = 2 ■ 1 19 atoms cm -2 ), exes: only for the 
surface column density: deviations where the expected surface 



TiRiFiC was tested by constructing a number of artificial ob- 
servations of a number of H i disks of a galaxy at a distance 
of 4 Mpc. Here we show the results of the fits, which were per- 
formed by an "observer" using TiRiFiC. In addition, we present 
the results derived with ROTCUR on idealised velocity fields, 
which were extracted from a artificial observation of the same 
disk with ten times lower noise and ten times higher velocity 
resolution. ROTCUR was run on a peak-velocity field and a 
first-moment velocity field, the input guesses reported in the 
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Fig. 9. Test 1 : Flat disk with constant scale height without noise. The orientation parameters, centre, systemic velocity and scale 
height were fitted "as one" with TiRiFiC. ROTCUR: input model with optimal guesses. Only the rotation velocity was fitted. For 
symbol- and colour coding see the introduction to this section. 
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column density is 0, i.e. where the artificial disk has already 
ended. 

(3) Method to determine the quantities, tir: use of TiRiFiC. ml : 
use of ROTCUR, fitting to the first-moment velocity field, pk: 
use of ROTCUR, fitting to the peak-velocity field. 

(4) Error in face-on Hi column density (10 18 atoms cm -2 ). 

(5) Ratio of error in face-on H i column density and expected 
face-on H i column density. 

(6) Error in rotation velocity (km s _1 ). 

(7) Error in inclination ( °). 

(8) Error in position angle ( °). 

(9) Error in right ascension of ring centre ( "). 

(10) Error in declination of ring centre ( "). 

(11) Error in systemic velocity (km s _1 ). 

(12) Error in scale height ( ")• 
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Table 1. Test 1: Flat disk with constant scale height without noise. The orientation parameters, centre, systemic velocity and 
scale height were fitted "as one" with TiRiFiC. Deviation in global dispersion: 0.01 km s _1 . ROTCUR: input model with optimal 
guesses. Only the rotation velocity was fitted. 
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Fig. 10. Test 2: Flat disk with constant large scale height; noise added to the cube. Orientation parameters, centre, systemic 
velocity and scale height fitted "as one" with TiRiFiC. In a first attempt the observer failed to achieve a good fit. The large 
scale height enhances drastically the ambiguity of inclination and rotation velocity. ROTCUR: input model with optimal guesses. 
Onlv the rotation velocity and the orientation carameters were fitted. For svmbol- and colour codine see the introduction to this 
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Table 2. Test 2: Flat disk with constant large scale height; noise added to the cube. Orientation parameters, centre, systemic 
velocity and scale height fitted "as one" with TiRiFiC. In a first attempt the observer failed to achieve a good fit. The large 
scale height enhances drastically the ambiguity of inclination and rotation velocity. Deviation in global dispersion: 0.57 kms -1 . 
ROTCUR: input model with optimal guesses. Only the rotation velocity and the orientation parameters were fitted. 
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Fig. 11. Test 2, after revision: same as 2, but this time started near to the correct result. The observer achieved a reasonable result, 
while TiRiFiC approached the x 1 minimum in parameter space very slowly. It is to be hoped that this problem will vanish after 
the implementation of a more effective x 2 minimiser. ROTCUR: input model with optimal guesses. Only the rotation velocity 
was fitted. For svmbol- and colour codins see the introduction to this section. 
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Table 3. Test 2, revised: same as 2, but this time started tear to the correct result. The observer achieved a reasonable result, 
while TiRiFiC approached the x 1 minimum in parameter space very slowly. It is to be hoped that this problem will vanish after 
the implementation of a more effective^ 2 minimiser. Deviation in global dispersion: 0.39 kms" 1 . ROTCUR: input model with 
optimal guesses. Only the rotation velocity was fitted. 



Quantity 


Range 


method 


A'hi 


A'hi 
Nai,exp 


Vrot 


i 


pa 


RA 


DEC 


Vsys 


Za 


u; 


(A) 




(A\ 

w 


(■>) 


(f\\ 
(O) 


\i> 


K°) 




(IV) 


fin 




A 


all 


tir 


160.22 


0.386 


4.0 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


a/a? 
V ZA 


all 
all 


fir 
111 


tj J.AL 


fl zL8^ 

U.40J 


u. / 


7 1 
/ . i 


u.u 


U. 1 


U. J 


n 1 

U. 1 


n 9 


Amax 


all 


tir 


1752.18 


0.975 


18.8 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 


70% 


tir 


8.43 


0.217 


1.1 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


a/a2 


/ U /C 


tir 
111 


1J.1Z 


U.il / 1 


1 .H- 


7 1 
/ . i 


u.u 


U. 1 


U.J 


fl 1 

U. 1 


fl 9 


A 

^max 


70% 


tir 


30.33 


0.608 


2.9 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 
ZA 


hi 
111 


tir 
Ill 




n ass 


1 1 

L.L 


7 1 




u.u 


U. 1 


U.J 


fl 1 

U. 1 


n i 

\j.L 


V i-i 


hi 


tir 
ill 


147 (S4 


S7^ 

U.J / J 




7 1 


u.u 


1 

U. 1 




1 

U. 1 




A 

^max 


hi 


tir 


361.90 


0.975 


6.6 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 


8-cr 


tir 


11.54 


0.478 


5.5 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 
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V ZA 


8 — IT 

o — u 


tir 
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1 8 S 1 

1 O.J 1 
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0. 1 
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/ . i 


u.u 
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U. 1 


U.J 
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U. 1 


fl 9 


A 


8 - cr 


tir 


48.88 
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18.8 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 


2 - <x 


tir 


6.70 


0.464 


6.3 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


a/a? 
V ZA 


9 — rr 


tir 
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lU.Ut 


u.jtu 


1 

y . l 


7 1 
/ . i 


u.u 
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U. 1 


U.J 


n 1 

U. 1 


n 9 
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2 - tr 


tir 


29.82 
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18.8 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 


0.5 - <x 


tir 


2.40 


0.523 


9.5 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


a/a? 
V ZA 


U. J C 


tir 
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^ S8 


U. JUiL 
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7 1 
/ . i 


u.u 


n i 

U. 1 


U.J 


n 1 

U. 1 


fl 9 


A 

ZA m ax 


0.5 - <t 


tir 


6.77 


0.720 


18.8 


7.1 


0.0 


0.1 


0.3 


0.1 


0.2 


A 


excs 


tir 


0.09 


















a/a? 
y ZA 


exes 


tir 
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tir 


0.18 


















A 


all 


nk 
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0.0 


0.0 


0.0 


0.0 


0.0 
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V ZA 


all 
all 
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u.u 


fi n 
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0.0 


0.0 
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0.0 


0.0 




A 
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oo 
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0.0 
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9 
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u.u 


u.u 


fl fl 
u.u 




Amax 
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P k 
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0.0 


0.0 


0.0 


0.0 


0.0 




A 


hi 


pk 






6.6 


0.0 


0.0 


0.0 


0.0 


0.0 






hi 


pk 






7.8 


0.0 


0.0 


0.0 


0.0 


0.0 




Amax 


hi 


pk 






12.0 


0.0 


0.0 


0.0 


0.0 


0.0 




A 


all 


ml 






9.7 


0.0 


0.0 


0.0 


0.0 


0.0 




VP 


all 


ml 






11.9 


0.0 


0.0 


0.0 


0.0 


0.0 




Amax 


all 


ml 






0.0 


0.0 


0.0 


0.0 


0.0 


0.0 




A 


70% 


ml 






5.9 


0.0 


0.0 


0.0 


0.0 


0.0 




VP 


70% 


ml 






7.6 


0.0 


0.0 


0.0 


0.0 


0.0 




Amax 


70% 


ml 






14.7 


0.0 


0.0 


0.0 


0.0 


0.0 




A 


hi 


ml 






13.0 


0.0 


0.0 


0.0 


0.0 


0.0 




VP 


hi 


ml 






14.3 


0.0 


0.0 


0.0 


0.0 


0.0 




A m ax 


hi 


ml 






20.1 


0.0 


0.0 


0.0 


0.0 


0.0 





G. I. G. Jozsa et al.: Kinematic modelling of disk galaxies, Online Material p 12 



- oo n oo rn oo s o 



-00 n 00 m 01 s 







0-^ 5r A # "0" # ig r 



50 200 250 



50 100 



200 250 



Fig. 12. Test 3: symmetric warp without central shift and without shift in systemic velocity. Noise was added to the data cube. 
ROTCUR: input model with optimal guesses. Only the rotation velocity and the orientation parameters were fitted. For symbol- 
and colour coding see the introduction to this section. 
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Table 4. Test 3: symmetric warp without central shift and without shift in systemic velocity. Noise was added to the data cube. 
Deviation in global dispersion: 0.03 kms -1 . ROTCUR: input model with optimal guesses. Only the rotation velocity and the 
orientation parameters were fitted. 
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/ 1 \ 


(2) 


(3) 


(4) 


/r \ 

(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


A 


all 


tir 


77.34 


0.505 


8.2 


5.4 


2.6 


0.3 


0.0 


0.1 


0.0 


VA Z 


all 


tir 


1/1/1 oz. 

144.86 


0.757 


12.5 


9.2 


1 a 

3.6 


0.3 


0.0 


0.1 


0.0 


A 


all 


tir 


471.18 


1.724 


32.1 


30.0 


9.5 


0.3 


0.0 


0.1 


0.0 


A 


70% 


tir 


18.61 


0.155 


2.5 


1.7 


1.1 


0.3 


0.0 


0.1 


0.0 
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T7 /if 
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0.201 


3.8 


2.2 


1 A 

1.4 
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0.3 


A A 

0.0 


A 1 

0.1 


A A 

0.0 


A 


70% 


tir 


58.22 


0.499 


7.9 


4.8 


2.5 


0.3 


0.0 


0.1 


0.0 


A 


hi 


tir 


64.09 


0.112 


2.1 


2.0 


3.3 


0.3 


0.0 


0.1 


0.0 


a/77 


hi 


tir 


85.23 


0.137 


3.2 


2.8 


4.5 


0.3 


0.0 


0.1 


0.0 


Amax 


hi 


tir 


204.22 


0.213 


7.3 


6.6 


9.5 


0.3 


0.0 


0.1 


0.0 


A 


8 - (T 


tir 


6.20 


0.815 


13.6 


10.3 


4.1 


0.3 


0.0 


0.1 


0.0 


V A 


C /T" 


tir 


1 AQ 
/ .60 


i 

I.Ujj 


1 A A 
10. 


i j. j 


A Q 


KJ.J 


u.u 


U. 1 


u.u 


^max 


8 - cr 


tir 


13.83 


1.724 


32.1 


30.0 


9.5 


0.3 


0.0 


0.1 


0.0 


A 


2 - (T 


tir 


5.32 


1.015 


16.4 


11.9 


2.8 


0.3 


0.0 


0.1 


0.0 


v A 2 


2 - cr 


tir 


6.77 


1.171 


18.7 


14.9 


3.1 


0.3 


0.0 


0.1 


0.0 


^max 


2 - cr 


tir 


12.93 


1.724 


32.1 


30.0 


5.1 


0.3 


0.0 


0.1 


0.0 


A 


0.5 - cr 


tir 


5.02 


1.149 


18.2 


13.2 


2.4 


0.3 


0.0 


0.1 


0.0 


V A 


U.j — cr 


tir 


A 

6.6 / 


1 OAO 
1 .ZuZ 


zu.U 


16. U 


O A 
Z.6 


U.J 


U.U 


n 1 

U. 1 


U.U 


^max 


0.5 - cr 


tir 


12.93 


1.724 


32.1 


30.0 


4.2 


0.3 


0.0 


0.1 


0.0 


A 


exes 


tir 


0.00 


















V A 


exes 


tir 


n nn 
u.uu 


















^max 


exes 


tir 


0.00 


















A 


all 


Dk 






16.0 


9.8 


1.4 


0.0 


0.0 


0.0 




\ A 2 


all 


pk 






21.5 


11.9 


2.0 


0.0 


0.0 


0.0 




A 


all 


Dk 






21.5 


20.7 


3.9 


0.0 


0.0 


0.0 




A 


70% 


nk 






6.0 


4.7 


0.2 


0.0 


0.0 


0.0 




a/72 

VA 2 


70% 


pk 






7.6 


6.2 


0.3 


0.0 


0.0 


0.0 




A 


70% 


Dk 






13.0 


10.2 


0.6 


0.0 


0.0 


0.0 




A 


hi 


pk 






18.2 


11.2 


1.0 


0.0 


0.0 


0.0 




VP 


hi 


pk 






23.1 


12.9 


1.4 


0.0 


0.0 


0.0 




Amax 


hi 


pk 






46.0 


20.7 


2.7 


0.0 


0.0 


0.0 




A 


all 


ml 






16.6 


12.4 


0.9 


0.0 


0.0 


0.0 




VP 


all 


ml 






19.9 


13.7 


1.5 


0.0 


0.0 


0.0 




Amax 


all 


ml 






37.7 


19.9 


3.5 


0.0 


0.0 


0.0 




A 


70% 


ml 






10.2 


9.0 


0.2 


0.0 


0.0 


0.0 




VP 


70% 


ml 






12.3 


10.1 


0.2 


0.0 


0.0 


0.0 




A 

^max 


70% 


ml 






17.5 


14.5 


0.3 


0.0 


0.0 


0.0 




A 


hi 


ml 






19.2 


12.9 


0.2 


0.0 


0.0 


0.0 




VP 


hi 


ml 






21.8 


13.6 


0.2 


0.0 


0.0 


0.0 




^max 


hi 


ml 






37.7 


18.9 


0.3 


0.0 


0.0 


0.0 
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Fig. 13. Test 4: asymmetric warp with shift of centre and systemic velocity; the scale height is kept constant with radius and noise 
was added to the data cube. ROTCUR: input model with optimal guesses. All parameters except expansion velocity were left 
variable. For symbol- and colour coding see the introduction to this section. 
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Table 5. Test 4: asymmetric warp with shift of centre and systemic velocity; the scale height is kept constant with radius and noise 
was added to the data cube. Deviation in global dispersion: below 0.01 km s _1 . ROTCUR: input model with optimal guesses. All 
parameters except expansion velocity were left variable. 



Quantity 


Range 


method 




A'Hi 

N Hl ,exp 


Vrot 


i 


pa 


KA 


DrlL. 


v sys 


Zo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


A 


all 


tir 


119.82 


0.126 


1.8 


1.2 


0.7 


0.9 


0.7 


0.5 


0.1 


VA 2 


all 


tir 


259.53 


0.248 


2.8 


1.9 


1.1 


1.5 


1.2 


1.0 


0.1 


A 


all 


tir 


965.85 


0.966 


8.7 


4.3 


3.4 


4.2 


4.3 


3.9 


0.1 


A 


70% 


tir 


17.63 


0.037 


0.6 


0.4 


0.3 


0.3 


0.2 


0.2 


0.1 


VA 2 


70% 


tir 


21.51 


0.046 


0.8 


0.4 


0.4 


0.4 


0.3 


0.2 


0.1 


A 


70% 


tir 


43.33 


0.088 


1.4 


0.8 


0.9 


0.9 


0.7 


0.5 


0.1 


A 


hi 


tir 


36.71 


0.083 


1.4 


0.8 


0.7 


1.0 


0.8 


0.6 


0.1 


v A 2 


hi 


tir 


56.58 


0.129 


2.1 


1.4 


1.1 


1.6 


1.3 


1.1 


0.1 


Amax 


hi 


tir 


187.26 


0.375 


5.4 


4.3 


3.4 


4.2 


4.3 


3.9 


0.1 


A 


8 - (T 


tir 


14.20 


0.192 


2.7 


2.0 


1.7 


2.1 


1.9 


1.4 


0.1 


v A 2 


8 - <r 


tir 


18.60 


0.223 


3.5 


2.5 


1.9 


2.5 


2.3 


1.9 


0.1 


Amax 


8 - cr 


tir 


36.53 


0.375 


5.4 


4.3 


3.4 


4.2 


4.3 


3.9 


0.1 


A 


exes 


tir 


6.45 


















VA 2 


exes 


tir 


6.45 


















Amax 


exes 


tir 


6.45 


















A 


all 


Dk 






34.0 


12.8 


0.4 


1.7 


1.1 


0.2 




V A 


oil 

all 


pk 






81. J 


lo.o 


V.y 


L. 1 


Z.U 


U.j 




Amax 


all 


Dk 






253.2 


39.9 


2.7 


in 


5.4 


1.4 




A 


70% 


Dk 






4.3 


7.2 


0.1 


0.6 


0.3 


0.1 




V A 


/u /c 








J.Z. 


O.U 


n 1 

U. 1 


u. / 




n 1 

U. 1 




A 

^max 


70% 


Dk 






8.6 


12.9 


0.2 


1.2 


0.6 


0.2 




A 


hi 


pk 






37.8 


12.7 


0.2 


1.0 


0.6 


0.1 




Va 2 " 


hi 


pk 






85.6 


16.9 


0.2 


1.3 


1.1 


0.1 




^max 


hi 


pk 






253.2 


39.9 


0.5 


2.9 


3.2 


0.3 




A 


all 


ml 






21.3 


11.7 


1.2 


2.8 


1.1 


0.8 




Va 2 " 


all 


ml 






34.9 


14.6 


3.5 


5.1 


1.8 


2.2 




^inax 


all 


ml 






97.9 


33.0 


11.0 


15.4 


5.3 


6.9 




A 


70% 


ml 






7.6 


6.7 


0.1 


0.9 


0.4 


0.1 




Va 2 " 


70% 


ml 






9.7 


6.9 


0.1 


1.0 


0.5 


0.1 




^max 


70% 


ml 






21.4 


9.2 


0.2 


1.5 


0.7 


0.2 




A 


hi 


ml 






20.5 


12.3 


0.1 


1.4 


0.6 


0.1 




VP 


hi 


ml 






35.5 


15.3 


0.2 


1.7 


0.7 


0.1 




A 

^max 


hi 


ml 






97.9 


33.0 


0.3 


3.2 


1.3 


0.3 
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Fig. 14. Test 5: asymmetric warp with central shift and shift in sytemic velocity. The scale height is constant with radius, the 
data cube contains no noise, constant scale height. The results get unreliable at low inclination (i < 25°). ROTCUR: first guess 
deviates by lOkms -1 in rotation velocity, by 10° in inclination and position angle, by lOkms -1 in systemic velocity, and by 10" 
in both central coordinates from the ontimal auess. All parameters exceot expansion velocity were left variable. For svmbol- and 
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Table 6. Test 5: asymmetric warp with central shift and shift in sytemic velocity. The scale height is constant with radius, the 
data cube contains no noise. The scale height is constant. Deviation in global dispersion: 0.06 kms 1 . The results get unreliable 
at low inclination (z < 25°). ROTCUR: first guess deviates by lOkms 1 in rotation velocity, by 10° in inclination and position 
angle, by lOkms -1 in systemic velocity, and by 10" in both central coordinates from the optimal guess. All parameters except 
expansion velocity were left variable. 
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22.9 


16.2 


36.1 


23.7 


16.7 




A 


70% 


pk 






13.9 


4.9 


0.5 


0.5 


1.7 


0.2 




Va? 


70% 


pk 






14.9 


5.8 


0.6 


0.8 


2.6 


0.4 




^max 


70% 


pk 






22.8 


9.8 


0.8 


1.9 


6.4 


0.9 




A 


hi 


pk 






22.7 


6.9 


0.7 


8.2 


4.7 


1.1 




Va? 


hi 


pk 






29.9 


8.4 


0.9 


15.1 


8.5 


2.2 




^max 


hi 


pk 






72.7 


17.1 


1.8 


36.1 


23.7 


6.1 




A 


all 


ml 






16.9 


7.3 


3.7 


7.7 


6.0 


3.6 




VP 


all 


ml 






20.0 


9.9 


9.6 


13.1 


9.8 


7.9 




Amax 


all 


ml 






42.6 


26.0 


30.1 


34.3 


22.2 


24.0 




A 


70% 


ml 






11.7 


4.2 


0.5 


1.9 


1.6 


0.5 




Va^ 


70% 


ml 






12.8 


4.9 


0.7 


3.2 


2.7 


0.9 




^max 


70% 


ml 






20.1 


7.0 


1.3 


6.7 


6.5 


2.2 




A 


hi 


ml 






14.1 


5.2 


0.8 


7.4 


4.5 


1.3 






hi 


ml 






15.6 


5.9 


1.0 


13.4 


8.1 


2.3 




Amax 


hi 


ml 






25.0 


8.7 


1.9 


34.3 


22.2 


6.0 





G. I. G. Jozsa et al.: Kinematic modelling of disk galaxies, Online Material p 18 




500 o cl o 




50 100 150 200 



50 100 150 200 



- 200 Q 



Fig. 15. Test 6: warp with central/sytemic velocity shift and variable scale height without noise added to the data cube. ROTCUR: 
first guess deviates by 5 km s _1 in rotation velocity, by 5° in inclination and position angle, by 5 km s _1 in systemic velocity, and 
by 5" in both central coordinates from the optimal guess. All parameters except expansion velocity were left variable. For 
svmbol- and colour codine see the introduction to this section. 
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Table 7. Test 6: warp with central/sytemic velocity shift and variable scale height without noise added to the data cube. Deviation 
in global dispersion: 0.06 km s _1 . ROTCUR: first guess deviates by 5 km s _1 in rotation velocity, by 5° in inclination and position 
angle, by 5kms 1 in systemic velocity, and by 5" in both central coordinates from the optimal guess. All parameters except 
expansion velocity were left variable. 
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Fig. 16. Test 7: asymmetric warp with central/sytemic velocity shift and variable scale height, noise added to the data cube. 
ROTCUR: optimal first guesses. All parameters except expansion velocity were left variable. For symbol- and colour coding see 
the introduction to this section. 
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Table 8. Test 7: asymmetric warp with central/sytemic velocity shift and variable scale height, noise added to the data cube. 
Deviation in global dispersion: below 0.01 kms" 1 . ROTCUR: optimal first guesses. All parameters except expansion velocity 
were left variable. 
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Fig. 17. Test 8: flat disk with constant scale height, centre, and systemic velocity. Noise was added to the data cube. ROTCUR: 
input model with optimal guesses. Only the rotation velocity and the orientation parameters were fitted. For symbol- and colour 
coding see the introduction to this section. 
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Table 9. Test 8: flat disk with constant scale height, centre, and systemic velocity. Noise was added to the data cube. Deviation 
in global dispersion: 0.12 kms -1 . ROTCUR: input model with optimal guesses. Only the rotation velocity and the orientation 
parameters were fitted. 
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Fig. 18. Test 9: slight symmetric warp without central/sytemic velocity shift and variable scale height, noise added. ROTCUR: 
input model with optimal guesses. Only the rotation velocity and the orientation parameters were fitted. For symbol- and colour 
coding see the introduction to this section. 
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Table 10. Test 9: slight symmetric warp without central/sytemic velocity shift and variable scale height, noise added. Deviation 
in global dispersion: 0.09 kms -1 . ROTCUR: input model with optimal guesses. Only the rotation velocity and the orientation 
parameters were fitted. 
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0.1 


0.0 


0.0 


25.7 


A 


70% 


tir 


12.01 


0.157 


1.5 


0.6 


0.4 


0.1 


0.0 


0.0 


1.5 


VA 2 


70% 


tir 


14.93 


0.232 


2.0 


0.7 


0.6 


0.1 


0.0 


0.0 


2.3 


A max 


70% 


tir 


33.38 


0.491 


3.7 


1.3 


1.4 


0.1 


0.0 


0.0 


5.7 


A 


hi 


tir 


32.72 


0.052 


1.3 


0.6 


0.6 


0.1 


0.0 


0.0 


0.6 


V A 


i->; 
ni 


tir 


ACl £,A 
4U.04 


n ri7^ 


1 s 
1 .o 


u. / 


U.o 


n i 

U. 1 


ft ft 
U.U 


U.U 


ft 7 
U. / 


A 


hi 


tir 


75.55 


0.183 


3.7 


1.3 


1.5 


0.1 


0.0 


0.0 


1.1 


A 


8 - (T 


tir 


9.78 


0.637 


6.2 


5.2 


2.7 


0.1 


0.0 


0.0 


5.2 




O — (T 


tir 


i 7 n7 

1Z.UZ 


n 7Q7 
U. ly 1 


7 ^ 


7 1 
/. 1 


J.V 


n i 
U.l 


U.U 


U.U 


O.U 


^max 


8 - cr 


tir 


25.56 


1.674 


16.2 


16.8 


9.7 


0.1 


0.0 


0.0 


9.2 


A 


2 - (T 


tir 


8.56 


0.786 


6.9 


6.4 


3.2 


0.1 


0.0 


0.0 


6.3 


V A 


Z — CT 


tir 


Q ^7 


n oni 
U.VU1 


Q ^ 
O.J 


o.U 


A A 
4.4 


n i 
U.l 


U.U 


n n 
U.U 


u.<5 


^max 


2 - cr 


tir 


16.74 


1.674 


16.2 


16.8 


9.7 


0.1 


0.0 


0.0 


9.2 


A 


excs 


tir 


8.25 


















A/A2 
V A 


exes 


tir 


o.jU 


















A 


exes 


tir 


10.00 


















A 


all 


Dk 






4.3 


8.5 


0.8 


0.0 


0.0 


0.0 




V A 


oil 

all 


pk 






A Q 


Q 7 


1 A 
1.0 


U.U 


U.U 


U.U 




A 

'-'max 


all 


Dk 






2.2 


15.0 


4.3 


0.0 


0.0 


0.0 




A 


70% 


Dk 






2.9 


5.6 


0.2 


0.0 


0.0 


0.0 




a/ A 2 
V A z 


/Uve 


pk 






j.j 


O. J 


n 7 
U.Z 


U.U 


U.U 


U.U 




A 

^max 


70% 


Dk 






5.7 


9.1 


0.3 


0.0 


0.0 


0.0 




A 


hi 


P k 






4.5 


7.7 


0.2 


0.0 


0.0 


0.0 




VP 


hi 


pk 






5.1 


8.9 


0.3 


0.0 


0.0 


0.0 




^max 


hi 


pk 






7.3 


13.7 


0.6 


0.0 


0.0 


0.0 




A 


all 


ml 






5.5 


10.3 


0.7 


0.0 


0.0 


0.0 




VP 


all 


ml 






7.6 


13.0 


1.4 


0.0 


0.0 


0.0 




Amax 


all 


ml 






19.0 


30.6 


3.8 


0.0 


0.0 


0.0 




A 


70% 


ml 






2.8 


6.4 


0.2 


0.0 


0.0 


0.0 




VP 


70% 


ml 






3.0 


6.6 


0.2 


0.0 


0.0 


0.0 




Amax 


70% 


ml 






4.6 


8.9 


0.3 


0.0 


0.0 


0.0 




A 


hi 


ml 






3.8 


7.9 


0.2 


0.0 


0.0 


0.0 




VP 


hi 


ml 






4.1 


8.1 


0.3 


0.0 


0.0 


0.0 




Amax 


hi 


ml 






5.7 


10.6 


0.5 


0.0 


0.0 


0.0 
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Fig. 19. Test 10: slight warp without central/sytemic velocity shift and variable scale height, noise added to the data cube. 
ROTCUR: input model with optimal guesses. Only rotation velocity, position angle, and inclination were fitted. For symbol- and 
colour coding see the introduction to this section. 
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Table 11. Test 10: slight warp without central/sytemic velocity shift and variable scale height, noise added to the data cube. 
Deviation in global dispersion: 0.12 kms" 1 . ROTCUR: input model with optimal guesses. Only rotation velocity, position angle, 
and inclination were fitted. 





Range 






N Hl ,exp 


^rot 




pa 


RA 


DEC 


^sys 


£o 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


A 


all 


tir 


89.66 


0.240 


5.6 


3.3 


1.3 


0.2 


0.0 


0.1 


3.6 


VP 


all 


tir 


215.03 


0.362 


13.3 


4.4 


1.6 


0.2 


0.0 


0.1 


5.1 


Amax 


all 


tir 


834.24 


1.000 


55.8 


9.7 


3.4 


0.2 


0.0 


0.1 


16.5 


A 


70% 


tir 


10.46 


0.094 


1.5 


1.6 


0.8 


0.2 


0.0 


0.1 


2.0 


V Zi 


IVj /C 


ur 


1Z.-4-D 


U. 1 ID 


1 Q 
1 .0 


1 Q 


\).y 


u.z 


u.u 


n 1 




^max 


70% 


tir 


25.22 


0.209 


3.2 


2.9 


1.6 


0.2 


0.0 


0.1 


4.0 


A 


hi 


tir 


21.19 


0.159 


2.5 


2.0 


1.3 


0.2 


0.0 


0.1 


2.4 


VP 


hi 


tir 


30.10 


0.229 


3.1 


2.4 


1.6 


0.2 


0.0 


0.1 


3.1 


^max 


hi 


tir 


65.39 


0.610 


6.7 


4.6 


3.0 


0.2 


0.0 


0.1 


7.0 


A 


8-0- 


tir 


11.92 


0.395 


8.3 


2.6 


1.6 


0.2 


0.0 


0.1 


3.8 


VP 


8 - cr 


tir 


16.09 


0.491 


18.0 


3.7 


1.9 


0.2 


0.0 


0.1 


4.3 


^max 


8-0- 


tir 


42.69 


1.000 


55.8 


9.7 


3.4 


0.2 


0.0 


0.1 


7.0 


A 


2-0- 


tir 


9.38 


0.520 


13.0 


3.7 


1.7 


0.2 


0.0 


0.1 


5.2 


VP 


2- o- 


tir 


10.33 


0.616 


25.0 


4.8 


2.0 


0.2 


0.0 


0.1 


5.3 




2-o- 


tir 


16.21 


1.000 


55.8 


9.7 


3.4 


0.2 


0.0 


0.1 


7.0 


A 


excs 


tir 


8.50 


















VP 


exes 


tir 


8.63 


















Amax 


exes 


tir 


10.00 


















A 


all 


pk 






4.1 


6.2 


0.4 


0.0 


0.0 


0.0 




VP 


all 


Dk 






6.3 


8.9 


0.7 


0.0 


0.0 


0.0 




Amax 


all 


P k 






16.0 


20.7 


2.1 


0.0 


0.0 


0.0 




A 


70% 


pk 






1.5 


2.7 


0.2 


0.0 


0.0 


0.0 




VP 


70% 


Dk 






2.3 


3.6 


0.2 


0.0 


0.0 


0.0 






70% 


P k 






5.5 


8.0 


0.3 


0.0 


0.0 


0.0 




A 


hi 


pk 






4.4 


6.8 


0.2 


0.0 


0.0 


0.0 




VP 


hi 


pk 






6.6 


9.3 


0.3 


0.0 


0.0 


0.0 




Amax 


hi 


P k 






16.0 


20.7 


0.6 


0.0 


0.0 


0.0 




A 


all 


ml 






3.2 


5.4 


0.8 


0.0 


0.0 


0.0 




VP 


all 


ml 






4.6 


7.2 


2.0 


0.0 


0.0 


0.0 






all 


ml 






6.1 


12.8 


6.4 


0.0 


0.0 


0.0 




A 


70% 


ml 






1.2 


2.6 


0.1 


0.0 


0.0 


0.0 




VP 


70% 


ml 






2.0 


3.7 


0.1 


0.0 


0.0 


0.0 




^max 


70% 


ml 






4.6 


8.8 


0.3 


0.0 


0.0 


0.0 




A 


hi 


ml 






2.4 


5.8 


0.2 


0.0 


0.0 


0.0 




VP 


hi 


ml 






3.6 


7.5 


0.3 


0.0 


0.0 


0.0 




A 


hi 


ml 






7.3 


12.8 


0.5 


0.0 


0.0 


0.0 
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Fig. 20. Test 11: symmetric warp without central/sytemic velocity shift, constant scale height. The data cube contains no noise. 
ROTCUR: input model with optimal guesses. Only the rotation velocity and the orientation parameters were fitted. For symbol- 
and colour coding see the introduction to this section. 
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Table 12. Test 11: symmetric warp without central/sytemic velocity shift, constant scale height. The data cube contains no noise. 
Deviation in global dispersion: 0.01 kms -1 . ROTCUR: input model with optimal guesses. Only the rotation velocity and the 
orientation parameters were fitted. 



Quantity 


Range 


method 


Mi i 


N Hl ,exp 


^rot 




pa 


D A 

KA 


DrlL 


v sys 


zo 


(1) 


(z) 


(1) 




f < » 
(•>) 


(6) 


(1) 


(8) 


(9) 


(1U) 


i 1 1 1 

(11) 


(1Z) 




all 


tir 


193.54 


0.629 


5.5 


4.4 


2.7 


0.0 


0.0 


0.0 


0.1 


A / a2 
V A 


oil 

all 


tir 


A?A 1 Q 

4j4. iy 


ft C7/1 
U.o /4 


i ft 1 
1U. 1 


A ft 
o.U 


A Q 
4. J 


ft ft 
U.U 


rt n 
U.U 


ft ft 
U.U 


ft 1 
U.l 


A 


all 


tir 


1489.84 


2.540 


32.6 


13.1 


11.6 


0.0 


0.0 


0.0 


0.1 


A 


70% 


tir 


9.06 


0.307 


1.6 


1.8 


0.7 


0.0 


0.0 


0.0 


0.1 


V A 


/Utc 


tir 


1 ^ QC 


ft 3"7C 
U.J /o 


Z.U 


Z.J 


i ft 

l.U 


ft ft 
U.U 


n n 
U.U 


ft ft 
U.U 


ft 1 
U.l 


A 


70% 


tir 


44.10 


0.685 


4.0 


5.0 


2.2 


0.0 


0.0 


0.0 


0.1 


A 


hi 


tir 


514.06 


0.497 


2.1 


2.4 


2.3 


0.0 


0.0 


0.0 


0.1 


VA 2 


hi 


tir 


719.46 


0.598 


2.6 


3.9 


3.7 


0.0 


0.0 


0.0 


0.1 


A 

u mas 


hi 


tir 


1489.84 


0.995 


4.8 


8.9 


9.4 


0.0 


0.0 


0.0 


0.1 


A 


8 - <r 


tir 


5.79 


0.448 


3.8 


4.2 


2.7 


0.0 


0.0 


0.0 


0.1 


-\l \2 
V A 


8 — cr 


tir 


1 J. J J 


ft ^"70 
U.j /Z 


A H 
4. / 


A Q 
4.^ 


/i ft 
4.U 


ft ft 
U.U 


rt n 
U.U 


ft ft 
U.U 


ft 1 
U. 1 


Amax 


8 - cr 


tir 


44.10 


0.995 


10.2 


9.2 


9.4 


0.0 


0.0 


0.0 


0.1 




2 - cr 


tir 


2.31 


0.370 


3.6 


3.9 


1.7 


0.0 


0.0 


0.0 


0.1 


V A 


Z — cr 


tir 


4.ZZ 


ft /i 

U.45Z 


A A 
4.0 


A H 
4. / 


Z.4 


ft ft 
U.U 


ft n 

U.U 


ft ft 
U.U 


ft 1 

U. 1 


^max 


2 - cr 


tir 


10.57 


0.902 


10.2 


9.2 


4.5 


0.0 


0.0 


0.0 


0.1 


A 


0.5 - cr 


tir 


1.48 


0.317 


2.5 


3.0 


1.3 


0.0 


0.0 


0.0 


0.1 


V A 


ft ^ ^i- 
U.J — (T 


tir 


o on 
z.yu 


ft /i 


Z.o 


1 A 
J. 4 


Z.U 


ft ft 
U.U 


ft ft 

U.U 


ft ft 
U.U 


ft 1 
U.l 


^max 


0.5 - <T 


tir 


8.12 


0.902 


4.0 


5.0 


4.5 


0.0 


0.0 


0.0 


0.1 


A 


excs 


tir 


0.02 


















V A 


exes 


ur 


ft fn 

U.Uj 


















A m „ 
"max. 


exes 


tir 


0.05 


















A 


all 


Dk 






17.3 


11.1 


1.0 


0.0 


0.0 


0.0 




V A 


all 

an 


pis 






ZJ.Z 


1 1 1 


i .j 


rt rt 

U.U 


ft ft 

U.U 


ft ft 
U.U 




A 

^max 


all 


Dk 






53.8 


25.4 


2.6 


0.0 


0.0 


0.0 




A 


70% 


nk 






5.4 


6.6 


0.3 


0.0 


0.0 


0.0 




V A 


I\)7C 


pk 






/ . j 


Q Q 
O.J 


ft /i 
U.4 


ft rt 
U.U 


ft ft 

U.U 


ft ft 
U.U 




Amax 


70% 


Dk 






14.6 


13.6 


0.6 


0.0 


0.0 


0.0 




A 


hi 


P k 






24.5 


11.0 


0.7 


0.0 


0.0 


0.0 




VP 


hi 


pk 






30.6 


12.7 


1.1 


0.0 


0.0 


0.0 




^max 


hi 


pk 






53.8 


19.7 


2.6 


0.0 


0.0 


0.0 




A 


all 


ml 






15.0 


8.7 


0.9 


0.0 


0.0 


0.0 




VP 


all 


ml 






23.8 


11.6 


1.2 


0.0 


0.0 


0.0 




Amax 


all 


ml 






22.4 


24.6 


2.1 


0.0 


0.0 


0.0 




A 


70% 


ml 






3.1 


4.0 


0.4 


0.0 


0.0 


0.0 




VP 


70% 


ml 






4.6 


5.3 


0.7 


0.0 


0.0 


0.0 




^max 


70% 


ml 






8.6 


8.9 


1.5 


0.0 


0.0 


0.0 




A 


hi 


ml 






22.4 


8.9 


0.8 


0.0 


0.0 


0.0 




VP 


hi 


ml 






29.2 


10.1 


1.1 


0.0 


0.0 


0.0 




A m ax 


hi 


ml 






53.5 


15.5 


2.1 


0.0 


0.0 


0.0 
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Fig. 21. Test 12: flat disk with constant scale height, noise added to the data cube. ROTCUR: input model with optimal guesses. 
Only the rotation velocity was fitted. For symbol- and colour coding see the introduction to this section. 
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Table 13. Test 12: flat disk with constant scale height, noise added to the data cube. Deviation in global dispersion: 0.08 kms -1 . 
ROTCUR: input model with optimal guesses. Only the rotation velocity was fitted. 



Quantity 


Range 


method 




%I 

N Hl ,exp 


Vrot 


i 


pa 


KA 




^sys 


zo 


/ 1 \ 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


/ 1 1 \ 

(11) 


(12) 


A 


all 


tir 


486.53 


0.678 


5.0 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


. / a 9 


all 


tir 


2020.45 


1.118 


7.6 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 


nil 
ail 


tir 

111 




? sin 


1 5 1 

i j . i 


n 4 


n n 

U.U 


n n 

u.u 


n n 

u.u 


n i 

U. 1 


n i 

U. 1 


A 


70% 


tir 


4.54 


0.198 


1.5 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


a/ A 2 


70% 


tir 


5.69 


0.325 


2.4 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 

n mas 


70% 


tir 
ill 


i n "so 




6 6 
u.u 


n 4 

U.t 


n n 

U.U 


n n 

u.u 


n n 

u.u 


n i 

U. 1 


n i 

U. 1 


A 


hi 


tir 


81.15 


0.086 


0.4 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A / A 9 


hi 


tir 


107. 39 


A 1 1 C 

0.1 15 


a 

0.6 


0.4 


A A 

0.0 


A A 

0.0 


A A 

0.0 


A 1 

0.1 


A 1 

0.1 


A 

^mas 


hi 


tir 


208.61 


0.242 


1.0 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 


8 - (T 


tir 


4.89 


0.951 


7.7 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A / A 9 

VA Z 


o 

8 - CT 


tir 


6.46 


1 TIC 

1.235 


9.2 


0.4 


A A 

0.0 


A A 

0.0 


A A 

0.0 


A 1 

0.1 


A 1 

0.1 


A 


8 - <T 


tir 


14.21 


2.810 


15.1 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 


2 - (T 


tir 


4.46 


1.052 


8.7 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


a/ A 9 

V A z 


Z — (T 


tir 


6.U4 


1.31 / 


a a 
9.9 


A A 
U.4 


A A 

u.u 


A A 

u.u 


A A 

u.u 


A 1 

U.l 


A 1 

U.l 


A 


2 - <T 


tir 


14.21 


2.810 


15.1 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 


0.5 - cr 


tir 


3.65 


1.062 


9.6 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A / A 9 

V A z 


U.D - (T 


tir 


4. / / 


1 a£/l 

1.364 


1U.0 


A A 
U.4 


A A 

u.u 


A A 

u.u 


A A 

u.u 


A 1 

U.l 


A 1 

U.l 


A 


0.5 - (T 


tir 


8.92 


2.810 


15.1 


0.4 


0.0 


0.0 


0.0 


0.1 


0.1 


A 


excs 


tir 


2.70 


















a / A 9 

VA Z 


exes 


tir 


4.06 


















A 

u aias 


exes 


ill 


X Q? 


















A 


all 








19.9 


15.2 


0.3 


0.0 


0.0 


0.0 




a/ A 2 

V A z 


all 


pk 






44.8 


23.7 


0.6 


0.0 


0.0 


0.0 




A 


all 


nk 






131.0 


56.5 


1.5 


0.0 


0.0 


0.0 




A 


70% 


nk 






1.9 


4.4 


0.0 


0.0 


0.0 


0.0 




a/ A 2 

VA Z 


70% 


pk 
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Fig. 22. Test 13: asymmetric warp with variable scale height, central shift and shift in systemic velocity. Noise added to the data 
cube. ROTCUR: optimal first guess. All parameters except expansion velocity were left variable. For symbol- and colour coding 
see the introduction to this section. 
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Table 14. Test 13: asymmetric warp with variable scale height, central shift and shift in systemic velocity. Noise added to the 
data cube. Deviation in global dispersion: 0.08 kms -1 . ROTCUR: optimal first guess. All parameters except expansion velocity 
were left variable. 
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Fig. 23. Test 14: asymmetric warp with variable scale height, central shift and shift in systemic velocity. Noise added to the data 
cube. The cloud flux used was 10~ 3 Jykms 1 for a radius below 48", 10~ 4 Jykms~' for a radius between 48" and 114", and 
1CT 4 Jy km s _1 beyond 1 14", correspinding to H i masses of 3.8 • 1O 4 M , 3.8 • 10 3 M o , and 3.8 • 1O 2 M , respectively, for a galaxy 
at a distance of 4Mr>c. the cloud masses scalins with the distance sauared. With this setuo. an inhomoeeneous distibution of 
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Table 15. Test 14: asymmetric warp with variable scale height, central shift and shift in systemic velocity. Noise added to the 
data cube. The cloud flux used was 1CT 3 Jykms 1 for a radius below 48", 10~ 4 Jykms -1 for a radius between 48" and 114", 
and 1(T 4 Jy km s _1 beyond 1 14", correspinding to H i masses of 3.8 • 10 4 M Q , 3.8 • 10 3 M , and 3.8 • 10 2 M , respectively, for a 
galaxy at a distance of 4 Mpc, the cloud masses scaling with the distance squared. With this setup, an inhomogeneous distibution 
of the H i was simulated. Since the orientation of the galactic disk is near to edge-on, no attempt was made to perform a fit to the 
velocity field. 



V^uaiiiiiy 


iVLUlgC 


IIlCLllOU 




JV H i,M'P 


Vrot 


i 


pel 


R A 


1 V 1 A 


v sys 


Zo 


en 


(7) 


\ J ) 




C51 


(6*) 


V 1 ) 




C91 








A 


all 


tir 


26.60 


0.471 


6.9 


8.2 


12.1 


7.9 


13.1 


3.7 


9.2 


VP 


all 


tir 


44.24 


1.013 


11.4 


10.0 


16.7 


13.2 


22.5 


4.7 


20.0 




all 


tir 


110.27 


4.126 


34.4 


17.5 


36.9 


35.9 


52.6 


10.2 


76.4 


A 


70% 


tir 


5.88 


0.130 


2.6 


5.0 


5.5 


2.1 


1.1 


2.1 


2.5 


•VA2 
V A 


l\) /c 


tir 


C ^1 
o.jZ 


n 1 £1 
U. 101 




u.U 


Q 1 

o. 1 


1 A 


1 C 
l.J 


z.o 


J.J 




70% 


tir 


21.57 


0.385 


4.2 


10.4 


18.0 


3.7 


4.0 


4.9 


6.6 


A 


hi 


tir 


42.58 


0.122 


2.8 


3.0 


2.2 


2.5 


0.7 


1.1 


2.6 


VP 


hi 


tir 


58.27 


0.142 


3.0 


3.5 


2.6 


3.6 


0.9 


1.3 


3.4 


Amax 


hi 


tir 


110.27 


0.252 


4.2 


6.4 


4.7 


8.4 


1.4 


2.4 


5.8 


A 


8-0- 


tir 


4.10 


0.795 


10.5 


9.1 


16.0 


12.9 


24.1 


4.6 


5.1 


VP 


8-cr 


tir 


5.76 


1.390 


15.4 


10.9 


20.6 


18.0 


31.1 


5.6 


5.7 


Amax 


8-0- 


tir 


11.99 


4.126 


34.4 


17.5 


36.9 


35.9 


52.6 


10.2 


8.2 


A 


2-o- 


tir 


3.79 


0.969 


12.1 


10.9 


19.0 


15.0 


29.9 


5.4 


5.1 


VP 


2-o- 


tir 


5.61 


1.554 


17.0 


12.1 


22.9 


19.9 


34.7 


6.2 


5.8 


Amax 


2-o- 


tir 


11.99 


4.126 


34.4 


17.5 


36.9 


35.9 


52.6 


10.2 


8.2 


A 


0.5 -cr 


tir 


2.36 


1.279 


17.7 


13.9 


28.0 


18.0 


41.1 


4.6 


7.1 


VP 


0.5 -a 


tir 


4.36 


1.933 


21.4 


14.4 


28.7 


21.9 


42.0 


5.2 


7.1 


Amax 


0.5 -o- 


tir 


11.99 


4.126 


34.4 


17.5 


36.9 


35.9 


52.6 


9.5 


8.2 


A 


excs 


tir 


0.00 


















VP 


exes 


tir 


0.00 


















Amax 


exes 


tir 


0.00 



















G. I. G. Jozsa et al.: Kinematic modelling of disk galaxies, Online Material p 36 




00 h 00 m 01?0 - 



42.5 
6D°OD'42'.'o 




-oo n oo m oi ro 



y - 0.2 



Fig. 24. Test 15: flat disk with constant scale height and solid-body rotation, no noise added to the data cube. While the fitting 
was in this case not performed without the knowledge of the parametrisation of the fake observation, convergence was reached 
very quickly: since TiRiFiC fits the surface brightness profile, the degeneracy of inclination and rotation velocity is broken. For 
svmbol- and colour codine see the introduction to this section. 
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Table 16. Test 15: flat disk with constant scale height and solid-body rotation, no noise added to the data cube. Deviation in 
global dispersion: 0.01 kms -1 . While the fitting was in this case not performed without the knowledge of the parametrisation of 
the fake observation, convergence was reached very quickly: since TiRiFiC fits the surface brightness profile, the degeneracy of 
inclination and rotation velocity is broken. 
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